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ABSTRACT 

As a step towards a better understanding of turbulent liquid-metal 
heat transfer in rectangular ducts with heat generation in the fluid 
stream, an analysis is performed for forced-convection heat transfer to 
slug flow in rectangular channels with aspect ratios from 1 to » and with 
heat sources in the fluid. The channel walls are uniformly heated, and 
the heat flux on the short sides is considered an arbitrary fraction or 
multiple of the heat flux on the broad sides. The analysis is based on 
the additional assumptions that (l) heat transport by eddy conduction is 
negligible compared with molecular conduction, (2) internal heat gener- 
ation is spatially uniform, and (3) fluid properties are invariant with 
temperature . 

The temperature distributions are determined by utilizing the method 
of superposition, and the required eigenvalues and constants are deter- 
mined analytically. The results obtained apply in the thermal entrance 
region of the channel as well as far downstream from the entrance. 

The effects of (l) the ratio which determines the relative role of 
internal heat generation to that of wall heat transfer, (2) specified 
aspect ratio, and (3) specified heat fluxes around the channel periphery 
on temperature distributions are investigated. Numerical results for 
wall temperatures and bulk-mean fluid temperatures are presented graph- 
ically. The solutions point out the locations of maximum temperature. 

The results are useful in estimating local heat-transfer characteris- 
tics in turbulent heat-generating liquid metal flow in rectangular ducts 
when the Prandtl and Reynolds moduli are low. 


INTRODUCTION 

Recent technological developments in space power generation have 
stimulated interest in the problem of forced-convection heat transfer to 
liquid metal flow in passages with internal heat generation in the fluid 
stream. This system has applications, for example, in the design of 
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liquid metal-fuel reactors, electromagnetic pumps and flowmeters, and 
liquid metal MHD generators. The fluid in these devices will be heated 
by radioactive fission products or by an electric current flowing through 
the fluid. A factor of importance for the proper operation of these de- 
vices is maintaining a satisfactory temperature distribution along the 
passage walls. The designer, therefore, must be able to compute the 
temperature distribution along the walls and to know how much heat must 
be removed to cool the walls to prevent temperatures from exceeding de- 
sign limits. The problem involves studying liquid-metal duct flow with 
combined internal heat sources and wall heat transfer. Turbulent flows 
are very often encountered in practice, and it is this flow regime which 
is of concern here. 

This investigation is concerned with hydrodynamic ally developed tur- 
bulent flow of a liquid metal in rectangular ducts with aspect ratios from 
1 to oo and with uniform internal heat generation. The channel walls are 
considered uniformly heated, but the heat flux on the short sides of the 
channel is an arbitrary fraction or multiple of the heat flux on the long 
sides. Such a heat transfer situation, for example, may be the result 
of unwanted heat leakage or addition through insulation. 

Attention is focused here on the rectangular duct because of its in- 
creasing use in- the applications mentioned. Within the knowledge of the 
author, the experimental and analytical studies of turbulent liquid-metal 
duct flows with internal heat generation have been confined to elementary 
geometries such as the circular tube (refs. 1 to 3, and 5) and the parallel- 
plate channel (refs. 4 and 5). In contrast to this moderate amount of 
information, turbulent-flow heat transfer to a heat- generating liquid metal 
in a rectangular duct has apparently received little analytical and no 
experimental work. A few studies related to the problem considered in the 
present investigation are noted. In the absence of internal heat genera- 
tion in the fluid, there have been developed solutions which approximate 
situations which might occur with liquid metals in turbulent flow through 
rectangular ducts. Fully developed slug-flow Nusselt numbers and wall- 
temperature distributions have been presented in references 6 to 8. In 
these references, the solution of the problem was obtained by assuming, 
in addition to a uniform velocity throughout the duct, that turbulent eddy- 
ing does not contribute to conduction of heat within the fluid. The results 
pertain to systems characterized by low Reynolds and Prandtl moduli and to 
the portion of rectangular duct beyond the thermal entrance region. In 
the discussion of reference 9, Hoagland discusses work done on the thermal 
entrance region for laminar slug flow in rectangular ducts, and reports 
some numerical results. Reference 10 has examined forced-convection heat 
transfer to laminar slug flow in a rectangular channel for the boundary 
condition of a duct wall temperature both peripherally and axially uniform. 
This analysis was carried out under the restriction of no internal heat 
sources. 
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It is the purpose of this investigation to study forced- convection 
heat transfer in a heat-generating liquid metal flowing in a rectangular 
channel where heating occurs on all four walls. The uniform heat flux on 
the short walls is assumed an arbitrary fraction or multiple of the flux 
on the broad walls. For the sake of completeness, the converse condition 
is also considered; namely, that the heat flux on the broad walls is any 
fraction or multiple of the flux on the short walls. Aspect ratios from 
1 to oo are considered for the rectangular channel, and the aim of the 
following analysis is the determination of the axial and peripheral 
temperature distribution and heat-transfer characteristics in the channel. 
The findings of the analysis should be applicable along the entire length 
of the ducts, that is, in thermal entrance as well as fully developed 
regions. 

The present very limited knowledge of turbulent liquid-metal flow and 
eddy-diffusivity-variation in noncircular passages make theoretical pro- 
gress very unlikely without appeal to simplified models. It does not seem 
likely, moreover, that any single model (for which a mathematical analysis 
is feasible) will prove adequate for all Reynolds and Prandtl moduli. 
Therefore, in order to gain some understanding of the complex problem of 
turbulent liquid-metal flow in rectangular ducts with wall heat transfer 
and internal heat sources, consideration will be given to a simplified but 
representative model; this specific model not only retains many of the 
physical characteristics of turbulent liquid metal duct flow, but also 
leads to a tractable mathematical problem. This model should, therefore, 
provide information on the temperature distribution and heat- transfer 
characteristics for such flows in rectangular passages. 

The idealized system assumed to approximate the forced convection 
system under consideration is based on the following postulates: (l) The 
established turbulent velocity profile is represented by a uniform 
distribution; (2) the thermal eddy diffusivity is small compared to the 
thermal molecular diffusivity and is neglected; (3) longitudinal heat 
conduction is small compared to longitudinal convection and transverse 
conduction and is neglected; and (4) the internal heat generation is spa- 
tially uniform. It is pointed out (e.g. , refs. 6 and 11) that the blunt- 
nosed turbulent velocity distribution for a liquid metal system can be 
represented satisfactorily by a uniform distribution. The second postulate 
implies (ref 12) that the thermal solution pertains to systems character- 
ised by low Prandtl moduli and low and intermediate Reynolds moduli. The 
third postulate has been shown in reference 13 to introduce a negligible 
error for Peclet moduli equal to, or greater than, approximately 100. 
References 6 to 8 point out that turbulent- flow heat transfer to liquid 
metals may be estimated, at least in the absence of internal heat gener- 
ation, by the use of the slug- flow solutions for molecular conduction. 

In this investigation numerical results are provided for the case of 
internal heat sources which are uniform across the duct cross section and 
along the duct length. The results can undoubtedly be extended, however, 
to include sources which vary in the transverse and longitudinal directions 
(refs. 14 and 15). 
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ANALYSIS 


The rectangular channel and its coordinate system are shown in fig- 
ure 1. . Turbulent velocity is postulated to he fully established at x = 0. 
The fluid temperature at the channel entrance is uniform across the sec- 
tion at a value t e > Within the channel a heating process takes place that 
includes a uniform heat generation within the liquid metal and a uniform 
heat transfer at the channel walls. The fluid is assumed to have constant 
physical properties, and only steady-state heat transfer is investigated. 

The established turbulent velocity profile is represented by 

u = U (1) 

The differential equation describing convective heat transfer for the 
idealized system takes the form 



The linearity of the energy equation (2) suggests that the temperature 
t(x,y, z) be written as the sum of two parts, 

t(*,y,z) = tq(x,y,z) + t q (x,y,z) (3) 

in which tQ corresponds to the situation where a heat-generating fluid 
with an entering temperature of zero flows through a channel with insulated 
walls (q B = q s = 0) and t q corresponds to the situation where a nongener- 
ating fluid entering at temperature t e flows through a channel with heat 
transfer q B and qg at the walls. The general solution can then be 
obtained by superposition of the two simpler solutions in accordance with 
equation (3) , since if the individual temperature fields satisfy the linear 
energy equation, then their sum does also. 


At this stage, it is convenient to employ dimensionless coordinates; 
the dimensionless equations and boundary conditions used to determine te 
and t q are then, respectively: ^ 


at Q = b 2 ^ + q^2 

^ By' 2 Bz' 2 * 

Bt Q 

=0 at y' = 0 and 1 (insulated walls) 


Bt Q 
Bz ' - 


(4a) 

(4b) 


at z’ = 0 


and a (insulated walls) 


(4c) 
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and 


■tQ(Ojy',z') = 0 (Entrance condition) 


5t q = + 

dx ' dy ' 2 dz ' 2 


( 4d) 


(5a) 


at q _ g B a _ 
ay' k 


2 k(o^+ a ) at ::y ' = 0 ( s P ec i f i e d wall heat flux) (5h) 


at 


q q_ 


ay’ 2 k(o + a) 


at y' = 1 (Specified wall heat flux) ' (5c) 


atq ^ q s a _ 
dz 1 


gq 


K 2k(o + a) 


at z 1 = 0 ( Specified wall heat flux) ( 5d) 


at. 


gq 


az 1 2k( o + g) 


at z 1 = cr ( Specified wall heat flux) ( 5e) 


■kq(0>ySz') = t e (Entrance condition) (5f) 

These two problems will be treated separately and the results combined to 
yield information for the general situation. There is first considered 
the problem of a heat- generating fluid flowing in an insulated channel. 

The problem of wall heat transfer to a nongenerating flowing fluid is then 
considered. 


Internal Heat Generation with Channel Insulated 

The solution for tq is found most easily by separate consideration 
of the fully developed and entrance regions. There is first considered the 
fully developed temperature which applies in the region downstream of the 
entrance region. The temperature tq^ satisfies equation (4a) 

St Q,d = *\,<1 a%,d q *2 

dx' Sy' 2 dz' 2 K 


( 6 ) 
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The fully developed situation for uniform internal heat generation is 
characterized "by the fact that 


_ Q . 

§x” PUCp ' ' 

Equation (7) states that the temperature at all positions in the rectangular 
cross section rises in the same linear fashion along the channel length. 

For the fully developed situation the boundary condition at the entrance 
of the channel (x =10) need not he considered, since it is accounted for 
by the entrance region solution, and equation (7) may be rephrased as 

5$”’ + < 8 > 

The function f(y',z') is found by inserting equation (8) into the differ- 
ential equation (6). This leads to the equation for f(y’,z') as 


d 2 f + d 2 f 
Sy ' 2 dz ' 2 


0 


(9) 


The boundary conditions on f(y',z') are determined from the thermal 
boundary conditions (eqs. (4b) and (4c)), so that 


at y' =0 and 1 (10a) 

at z 1 = 0 and a (lOb) 

Taking a solution of equation (9) in the form 

f(ySz') = aQy’2 + a jy' + a 2 z f2 + a^z T + a 4 (ll) 

in which aQ, a-]_, ... , are constants chosen to satisfy the boundary condi- 
tions (eqs. (10a) and (10b)), there is obtained the result 

f(yS z T ) = a 4 “ constant (l2 a ) 

The constant can be evaluated from an overall energy balance on the fluid 
for the length of channel from 0 to x : 


bf 


= 0 


df 

’5F’ 


= 0 
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This leads to the condition from equation (8) : 



fCy’jZ^dy'dz' 


G 


(12b) 


from which the constant is evaluated. The resulting expression for 
f(y' , z') is found to he 


f = 0 


(13) 


The final expression for. the temperature distribution t^^ is given by 


JSd _ x . 

Qa 2 /k ~ 


(14) 


which applies only in the fully developed region downstream of the thermal 
entrance region. 


To determine the temperatures in the entrance region, it is conven- 
ient to introduce an entrance temperature tq so that 

t Q = t Q,d + bQ ( 15 ) 

From the linearity of the energy equation (4a), it is found that to 
satisfies the equation 


dt$ S 2 tQ S 2 tQ 
Sz 1 Sz T ^ 


with the boundary conditions 

St Q = 

V 

OtQ ^ 

dz ’ - 


0 


0 


at y' = 0 and 1 


at z T = 0 and a 


(16a) 


(16b) 


(16c) 
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At x* = 0, the condition is 

tQ(0,y',z') = 0 = tQ^ d (0,y T ,z’) + t^(0,y',z’) 
or, hy rearranging: 


tQ(0,y',z') 

Qa 2 //c 


= -f(y'#z') = o 


(led) 


The solution of equation (16a) which will satisfy equations (16b) to 
(I6d) can he found hy using a product solution which leads to a separation 
of variables. This will have the form 


00 00 


t 


Q 


Qa 2 /k 


22 - 


-A? -x’ 
Y(y')Z(z')e 1J 


(17) 


i=0 j=0 

in which the functions Y(y') and Z(z') satisfy the differential equations 


d 2 Y 

dy’ 2 

d 2 Z 


+ cx?Y = 

- & - 


0 


0 


(18a) 

(18b) 


with the respective boundary conditions 

at y' = 0 and 1 (18c) 


at z' = 0 and a (l8d) 

2 2 2 

and A d j = cc^ + J3j, Equations (18 a) and (18b) with their boundary condi- 
tions belong to the well-known class of differential equations of the 
Sturm- Liouville type. Solutions are possible only for discrete, though 
infinite, sets of cy_ and Pj values. 

It may be seen readily that the differential equations (18a) and (18b) 
and the boundary conditions are satisfied if there are taken for the func- 
tions Y(y') and Z(z’) the expressions 




dY 

dy' 

dZ 
dz ' 


= 0 


= 0 


Y = cos(iity ') ; 



i = 0, 1, . . . 
j = 0, 1, . . . 


(19 a) 
(19b) 
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The eigenvalue Aj_j is then given by the following expression: 




4 2 + i 


(19 c) 


The coefficients a-jj in equation (17) are evaluated to satisfy the 
entrance boundaiy condition equation (l6d). This gives the condition 


CO 00 




j*z’\ t Q (o,y',z’) 


Qa 2 /* 


i=0 J=0 


= 0 ( 20a) 

It follows immediately that 

a ij = 0 for a11 1 and J (20b) 

This result indicates that there is no entrance region solution; that is, 
t^(Qa 2 //c) =0. 

Now that t<^ d and t^ are known, they can be superposed as in 

equation (15) to obtain the solution that applies over the entire length 
of the channel, which is simply w 


Qa 2 ' /k 


= x* 


( 21 ) 


The local bulk fluid temperature tQ ^(x') along the channel length, for 
the uniform heat source, is given by ’ 



x 


or 





( 22 ) 


Wall Heat Transfer Without Internal Sources 

There is considered next the situation where there are prescribed 
wall heat fluxes q B and q s = aqp at the channel walls (fig. 1) but no 
internal heat sources. The temperature t q (x»,y*, z *) is the solution to 
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equation (5) . To obtain a solution fox* t q that will apply over the entire 
length of the channel, it is convenient to break t q into two parts. The 
first part is t q ^ 3 , the fully developed solution. The second part is tq, 
an entrance region solution that is added to vtq 3 to obtain temperatures 

in the region near the entrance of the channel. The temperature tq is 
then given by 4 

= ^q, d + ^q 

From the linearity of the energy equation, t q 3 and t q have to each 
satisfy equations ( 5 ). 9 

The fully developed solution is considered first. Far from the en- 
trance of the rectangular channel the temperature rises linearly in the 
axial direction because of the uniform, but unequal, heat inputs at the 
channel walls. From a heat balance on the fluid, the temperature gradient 
in the fully developed region must be 


4 = constant (24) 

pUcpa^a 


The temperature distribution is then given by the equation 



q A 




+ F(y ' 



(25) 


An equation for the function F(y'>z') is found by substituting equation (25) 
into equation (5a). This gives 


d 2 F + S 2 F 
by ' 2 c!z ' 2 


(26a) 


The boundary conditions on F(y , ,z I ) are determined from the thermal boundary 
conditions (eqs. (5b) to (5e)) so that 


3F 

dy' 


2(a + aj 


at y' = 0 


Sf 

"Sy 7 

Bf 
bz f 

SF 
6z ' 


2(a + a) 


at 


aa 


2 (a + a) 


at 


acr 


2(a + a) 


at 


y’ = i 


z' = 0 


Z' = CT 


} 


(26b) 


(26c) 
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The solution of equation (25a) is taken in the form 

F(y'jz') = t>oy' 2 + *17' + ^2 Z ' 2 + + (27) 

The constants Tdq, b B , hg, and hj can he found from the conditions given 
by equations (26b) and (26c) and which yield 

*0 - afaV V T 131 = ' t ’°' ^ = 2(* + W b5 ' ' c,1>2 

Hence equation (27) can be represented in the following form: 

F(y T ,z') = 2 (a ' + a) [j' 2 “ y ' + a Z ’ 2 " aZ _D + '• ( 28a ) 

The remaining constant b^ can he evaluated from an overall energy balance 
on the fluid for the length of channel from 0 to x : 



t e ) dy 1 dz 1 


q 

pUcpa 2 


X 


This leads to the condition from equation (25): 



zOdy'dz’ = 0 


(28b) 


from which the constant is evaluated. The expression is then inserted 
into equation (25) to give the final expression for the fully developed 
temperature distribution: 


Sul I S - i v> 

q/K " a 


+ 


1 

2(a + a) 



- y’ 




(29 a) 


This can be rephrased in an equivalent form 


~kq, d ~ ^e 
q B a//c 



+ y ' 2 



az 




a< 


a) 


Equation (29h) applies only in the fully developed region. 


(29b) 


To determine the temperatures in the thermal entrance region the 

function t* is needed. From the linearity of the energy equation, the 

^ * / \ 
equation for t is the same as equation (5a): 

Si 

at* a 2 t* i a 2 t* 

dx ' dy ' 2 d z ' 2 


(30a) 
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Since the wall heat additions have already "been accounted for in the fully 
developed solution, the .boundary conditions for tq are that no heat is 
transferred at the walls: 4 


atq 


= 0 at y' =0 and 1 


K 


■^T = 0 at z ' =0 and a 

At the channel entrance (x T = 0) the condition is 

tq(0,y',z') = t e = t q ^ d (0,y',z') + t*(0,y',z') 


(30b) 

(30c) 


or by rearranging. 


tq(0,y',z') = - _ 


, , a ,2 
y' + — z - az 


1 + | (l + aa)j 


( 30d) 


It will now be convenient to represent the entrance region temperature tq 
of the fluid by two functions e(x',y') and ^x^z 1 ) such that y 




= e(x',y') + |(x',z') 


In terms of 6 the energy equation (30a) becomes 

Se _ 5 2 e 

By'2 

for which the boundary conditions will be taken as 

Be 

■^rr = 0 at y' = 0 and 1 


0(0, y’) = - y*2 - y* + 1 (i + aa)j 


Correspondingly , the function \[r is given hy. 


d\J f _ 




with the boundary conditions 




Sz 1 


= 0 at z T = 0 and a 

*(0,z') = - ^ z* 2 - az*j 


(31) 


(32) 


(33) 
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A solution for the function 0 is taken in the form 


0 = 


00 

A m e _rmX cos (r^y’) 


m=2, 4, . . . 

in which the eigenvalues r m are given by 


( 34 ) 


r m = nw 


In a similar manner the function t has the solution 


(35) 


♦ 



cos (SjjZ 1 ) 


where the eigenvalues 


6 n are given by 


(36) 


&n - ™ (37) 

The coefficients A m of equation (34) and the coefficients Bn of 
equation (36) are evaluated to satisfy the respective entrance boundary 
conditions. This gives the conditions 



Am cos (mny ’ ) = - £y ' 2 . y . + | (l + a af] 


(38a) 


00 

I 


B n cos 


- az 


■) 


(38b) 


( n *v)-- (f z ’ : 

n=2,4, ... 

According to Sturm-Liouville theory the coefficients A m and ^ given by 


% = - 


^y'2-yi + — (l + aa ^J cos (nijfy * ) dy 1 


£ 


cos 2 (mny’Jdy' 


(39 a) 
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( 391 :) 


The integral appearing in the numerator of equation (39a) may he written as 



i (l + au) cos (mity^dy' 

6 — i 


nr jt* 


where m is an even integer. If m is an odd number, the integral yields 
a value of zero. The series coefficients are thus 


A-m = - 


2 2 
m jt 


} m = 2, 4, 


(40a) 


In a similar manner the coefficients B n are obtained as 

4clc7 

^ = - -s-s* n = 2, 4, . . . (40b) 

n^jT 

Now that t„ j and t* are known, they can be superposed as in 
q 1 

equation (23) to obtain the solution which applies over the entire length 
of the channel. This is 


— 7—^ sc 2(1 + + y T ^ - y ? + — z*2 . az T + i (l + aa) 

q B a/K \ a) ^ J o 6 


00 

z 


— cos (roity ' ) e” m x 


m=2, 4, . , . 


a a > - cos 

n^jt^ 


( z ’ \ 

[mt — I 

{ °) 


zj_\ -(n 2 jt 2 x'/a 2 ) (41) 


The local bulk mean temperature tq^ along the channel is given. by 


15 


: q,l> 


= t 0 + 


ptfd p aa 2 


X 


= t e + — x' 
e OK 


Substituting q = 2 (a + a)(q p a) into equation (42a) gives 


t q,t. - t e _ 

q B aA " 



(42a) 


(42b) 


The analysis of the preceding paragraphs led to the determination of 
the temperature distribution for the condition where the heat flux on the 
short sides is an arbitrary fraction or multiple of the heat flux on the 
broad sides, or more simply, where q g = aq B . There may arise, however, 
practical applications where it is desirable to express the heat flux on 
the broad walls as q B =. aq g . The temperature distribution in this situation 

is then obtained from equation (41) by substituting a = 1/a and q B = aq s : 



m=2, . • 


y') + ^ z '2 - s' + i (a + o) 

-fg cos 


- a 


00 

2 ss c ° s (”" 


/a2 


n=2, , 


The local bulk mean temperature is given by 

t q^t " t e J- i\ . 

q s a//c \ a) 


(43) 


(44) 


Combined Internal Heat Generation and Wall Heat Transfer 

Results for the situation where internal heat generation and wall heat 
transfer are occurring simultaneously are f$und by combining the solutions 
for t Q and t q in accordance with equation (3). For side wall heat 
transfer given by q s = aq B , the solution can be “written as 
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where 


R = — 


% 


(46) 


the temperature distribution may also he expressed in terms of side wall heat 
transfer by substituting l:/ci. for a, ecq g for q-g, and replacing the 
parameter R by the parameter S, defined as 


S 


Qa 


(47) 


The parameters R and S are the ratios of internal heat evolution to the 
heat transferred at the channel walls and give a measure of the relative 
importance, in connection with temperature development, of internal heat 
generation in the presence of wall heat transfer. 


Heat Transfer Results 

From the temperature distribution given by equation (45), various 
quantities of engineering interest can be determined. In the developments 
to follow, the analytical results will be expressed in terms of the broad 
wall heat transfer q-g. 

Results of practical interest are the wall temperature variations 
corresponding to prescribed wall heat transfer and internal heat genera- 
tion. The local temperatures t v (x T ,z ! ) = t v ^ along the broad walls can 
be found from equation (45) by evaluating it at y T = 0 or at y* = 1: 
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n= 2 , . . . 


The local temperatures t w (x',y') = t ws along the side walls are found in 
a similar manner from equation (45) by evaluating it at z' = 0 or at 
z' = at 


^ws ~ 



+ y ' 2 - y' + ^ (1 + aa) 



4 , -m 2 jr 2 x' 

— 5-9 cos (mity')e 


aa 


00 

1 


-(n 


2 2 


7° 2 ) 


2 2 


(49) 


n= 2 , < 


Another form of these equations which is more convenient is obtained 
by introducing the bulk mean temperature For a uniform heat source 
and uniform wall heat transfer, the bulk temperature is given by 


*b - *e 



(50) 


Then the local wall- to bulk-temperature difference along the braod walls 
is given by 
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c wB t- b a ,2 , , 1 / 

- z - ccz' + - (1 + aa) - 


f/c a 


4 -m 2 it 2 x T 
2 2 e 


( n * # 


(n 2 n 2 x'/cr 2 ) 


n — 2, . . . 


while the temperature difference along the side walls is 


^ws • ^b .9 t , 1 
— * — = y' c ~ y 1 + -=• 


(l + aa) - 


2 2 


(mity 1 ) 6 




Tn — 2, . . 


4 - (n^n^x 1 /a^) 

2 2 e 


n=2, . . . 


A noteworthy feature of equations (51) and (52) is that the variation of the 
local wall- to hulk- temperature differences is independent of the internal 
heat generation rate. The differences between the fully developed wall and 
hulk temperatures are 


(^wB - ^h) 

% a A 


a , 2 T , 1 

— z 1 - az T + — 

0 6 


T + i (l + do) 

6 


(53 a) 


(t W s - 't'b)(3 o 1 , . 

— i K = y ' 2 - y’ + s (i + ao) 


(53h) 


The ratios of local to fully developed temperature differences along 
the broad walls or side walls at any location in the channel is found from 
equations (51) to (53) as 


t IwD- ’t-b) _ m=2, ■ 

(^wB _ ^b) (j 


4 -m 2 rt 2 x’ \ 4 

m 2 n 2 e +(W ) n 2«2 COS 


n=2, . . . 


( njr v) ! 


- (n 2 rt^x f /cr 2 ) 


— z ! ^ - az f +•£■(! + a a) 


(54a) 
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(t-ws " ^b) 

( "t-ws _ tb ) ^ 



cos (uuty'je 


y' 2 - y' + i (l + acr) 


+ aa 


oo 

I 


n=2 , . . . 


2*2 


n^jt 


-(n 2 jt 2 x T /a 2 ) 

e 


(54b) 


It may be observed from the foregoing results that a prescribed heat flow 
from the duct surfaces to the fluid produces, at a given axial location, a 
local wall temperature or wall- to bulk- temperature difference which 
varies around the periphery of the duct and which also varies with duct 
aspect ratio. 


Results and Discussion 

To illustrate the effects of internal heat generation and prescribed 
wall heat transfer on wall temperature distributions, local bulk mean 
temperature, and thermal entry lengths for liquid metal flow in rectangu- 
lar passages, a number of solutions have been obtained for various 
combinations of aspect ratio a, wall heat flux ratio a or a, and heat 
flux ratio R or S. For the sake of brevity, only some of these results 
are included in this paper. In particular, values of the parameters a 
and a chosen for the computations correspond to the following cases: 

(a) a = -1, in which the broad walls are heated and the side walls are 
cooled; (b) a = 0, in which the .side walls are insulated; (c) a = 1, in 
which uniform heating (or cooling) takes place all around the duct periph- 
ery; and (d) a = 0, in which the broad walls are insulated. 

Designers of channels for the applications previously mentioned are 
interested in the temperatures achieved by the walls, and in particular 
the peripheral location where the wall temperature will assume its high- 
est value for a known wall heat input and internal heat generation rate. 
From an examination of the analytical results, it is to be expected that 
the peak temperatures will occur either in the corner of a duct or else 
at the centerline of the broad wall or short wall, depending upon the wall 
heat flux ratio a. Therefore, knowledge of temperature conditions in 
such regions is of special interest. The wall temperatures are given 
here relative to the bulk temperature, since the heat flux ratio R or 
S is then eliminated as a parameter. 

The longitudinal variation of the dimensionless wall temperatures 
were evaluated from the analytical solution and are presented in fig- 
ures 2 for duct aspect ratios of 1, 4, and 10. For the special case of 
side walls insulated (a = 0) the solution is independent of the duct 
aspect ratio, and therefore, the result applies for all aspect ratios. 

At some places along a wall, the wall- to bulk-temperature difference 
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■may be negative. This is understandable, however, if it is recalled that 
W is a local value along the wall, while t b is an average value over 
t e entire cross section. When the broad walls are heated and the side 
walls are cooled, the peak temperatures occur along the centerline of the 
broad walls at all axial positions for the aspect ratios shown (fig. 2(a)) 
With heating from only the broad walls (fig. 2(b)), these walls assume a 
uniform temperature for all aspect ratios which is higher than the tempera- 
tures along the insulated side walls. For a rectangular duct with uniform 
heat flux all around the periphery (fig. 2(c)), the peak temperatures for 
all a s occur at the duct corners. Finally, when heat is transferred 
from only the short walls (fig. 2(d)), these walls attain a uniform temper- 
ature which is higher than the temperatures along the adiabatic broad walls. 

Another quanitiy which is of practical interest is the bulk mean tem- 
perature variation along the length of the channel. The bulk temperature 
is given here relative to the temperature of the fluid at the entrance to 
the channel. The dimensionless bulk temperatures are presented in figures 3 
for duct aspect ratios of 1, 20, and °° (parallel-plate channel) with the 
heat flux ratio appearing as a family parameter. 

Positive and negative values of the parameters R and S are con- 
sidered in the figures. It is supposed that Q is positive (a heat source). 

Ju P ? S J t:i T e ValUe ° f R ' therefore > implies that q B is positive, that is, 
that heat is being transferred from the broad walls to the fluid. A 
negative value of R, on the other hand, implies that q B is negative or 
that heat is being transferred from the fluid to the broad walls. For 
positive R, therefore, internal heat generation and broad-wall heat transfer 
reinforce one another to produce a bulk temperature larger than that obtained 
in the absence of internal heat generation. Conversely, for negative R, 
the broad wall heat transfer opposes internal heat generation in the bulk 
temperature development. Similar arguments apply in connection with side 
wall heat transfer q s and heat flux ratio S. 

v° r vei ^ r STn all values of |R| or J S| , wall heat transfer dominates 
the bulk temperature development, while for large values of |R| or 
I Sj , the effects of internal heat generation dominate. This accounts for 
the varieth of trends that are evident in each of the figures. It is also 
observed that the bulk temperature development is, in general, slightly 
affected by duct aspect ratio cr to about 20 and insignificantly thereafter. 
For a duct with insulated side walls, however, the bulk temperature de- 
velopment is independent of the duct aspect ratio. 

The foregoing presentation of results has been concerned with wall 
temperatures in the thermal entrance region. As a matter of general 
interest the wall temperature variation around the periphery of the duct 
in the fully developed region is considered. The wall temperatures are 
again given relative to the bulk temperature, since in the fully developed 
region the temperature difference (t w - tb)^ is independent of x. In 

addition, 'as noted earlier, the temperature difference t w - tb is inde- 
pendent of the heat flux ratio R at all axial positions. 
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Fully developed vail temperatures are presented in figures 4 for 
values of a and a considered earlier and for various aspect ratios. 

It is worth while to recall that t w is a local value along the wall, 
while t*b is an average value over the entire cross section. Therefore, 
at some places along the wall, the temperature difference (tv - tb)^ is 
negative, which means that tb is larger than tv The hot spots are 
strikingly displayed in these figures, appearing in the corners, or at the 
broad wall or short wall midpoints. It is also seen that the aspect ratio 
has a profound effect on the temperature distribution. 

Of considerable practical importance to the designer is the knowledge 
of the conditions under which entrance effects must be accounted for in 
heat-transfer calculations. The approach of a local wall- to bulk- 
temperature difference to the fully developed value is, in theory, asymp- 
totic. Consequently, it is difficult to identify a specific length of 
channel as a thermal entrance length. It is practice to define a thermal 
entrance length in terms of the downstream distance x/aRePr at which 
the temperature difference approaches to within 5 percent of the fully 
developed value. The variation with dimensionless axial distance of local 
to fully developed wall temperatures at the duct corner and at the broad 
wall and short wall centerlines has been evaluated from the analytical 
expressions. The results thus obtained have been plotted in figures 5 for 
values of a = - 1 and u = 0, for parametric aspect ratio values. It 
should be noted that the corner temperature ratio is not shown for the 
square duct when a = - 1, since for this situation the corner temperature 
difference is zero at all axial positions. Lines delineating the condition 
(t v - t*b)/(tv - tt,) d = 0.95 have been drawn In the figures to facilitate 
determining the thermal entrance length. 

From an Inspection of the graphs, it is seen that for the wall heating 
conditions represented, the wall temperature profiles become fully de- 
veloped at widely different distances from the entrance to the heating sec- 
tion, depending upon the particular wall location chosen for consideration. 
It is also evident that there Is a strong influence of the side walls of 
the ducts on the temperature developments, with the thermal entrance 
lengths Increased with increased aspect ratio. 

It is practice in reporting heat transfer connected with flow through 
noncircular passages to present average heat transfer coefficients or 
Nusselt numbers based upon a heat flow averaged around the duct periphery 
and on an average wall temperature. This practice has utility when the 
wall temperature remains constant everywhere, or is at least constant 
around the periphery of a duct at a given axial position. In the present 
situation where the wall boundary condition is one of peripherally and 
axially uniform heat input, however, it is apparent that the knowledge 
of the resulting local wall- temperature distribution Is of more Importance 
to the designer than the Nusselt numbers or even average thermal entry 
lengths. Therefore, no attempt has been made to determine these quantities. 
It is felt that the results are presented in a form more convenient for 
engineering calculations. 
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CONCLUDING REMARKS 

Results nave been presented for heat transfer in rectangular channels 
with prescribed wall heat fluxes and heat sources uniformly distributed in 
the fluid. The various effects considered have thrown some light on the 
axial and peripheral temperature distributions as well as on the heat 
transfer characteristics. 

The present results are strictly valid only for the slug- flow velocity 
distribution with heat transfer only by molecular conduction. The simpli- 
fication provided by these assumptions, however, has made it possible to 
obtain exact mathematical solutions to the governing energy equation. Of 
greater importance is the fact that these temperature distributions approxi- 
mate conditions to be expected for turbulent liquid-metal flow in rectangu- 
lar ducts for relatively low Prandtl and Reynolds moduli. The results, in 
addition, point out the locations of maximum temperatures. 

For improved heat transfer calculations, velocity-profile and eddy 
diffusivity variations would have to be taken into account to provide a 
more realistic description. . Turbulent velocity profiles for circular pipe 
and parallel plate duct systems have been satisfactorily represented by 
power-law expressions. Within the knowledge of the author, velocity-profile 
and eddy diffusivity distributions for turbulent flow in rectangular ducts 
have received little theoretical consideration. Reference 16 presents a 
variational method for determining velocity distributions for flow of a 
power law fluid in cylindrical ducts. The results might prove useful for 
heat transfer studies. 

In closing, it should be mentioned that the slug-flow, molecular- 
conduction analysis can be used to treat other wall, boundary conditions, 
such as a specified axial and peripheral wall temperature distribution. 

NOMENCLATURE 

Am coefficient in series for'; temperature distribution with wall 

heat transfer and Q = 0 

length Of short side 

coefficient in series for temperature distribution with internal 
heat generation and qp = Is = 0 

coefficient in series for temperature distribution with wall 
heat transfer and Q = 0 

b length of broad side 

. . ’ 

Cp specific heat of fluid at constant pressure 


a 

a ij 

Bn 


F(y',z') 

f(y';Z') 

Pr 

Q 

q 

q 

R 

Re 

S 

t 

U 

u 

X 

x’ 

Y(y’) 

y 

y’ 

Z(z») 

z 

z* 
z« * 

a 


transverse temperature distribution in fully developed region 
with wall heat transfer and Q = 0 

transverse temperature distribution in fully developed region 
with internal heat generation and qg = q s = 0 

Prandtl number, M-Cp//c 

rate of internal heat generation per unit volume 
heat addition per unit wall area 

heat addition per unit channel length, equal to 2bqg + 2aq s 
heat flux ratio, Qa/q^ 

Reynolds number, pUa/p 
heat flux ratio, Qa/q s 
temperature 

velocity taken as uniform over duct cross-section 
local fluid velocity 

coordinate measured along the axial direction 
dimensionless axial coordinate, x/aRePr 

eigenfunction for case of internal heat generation and q^ = q g 
coordinate measured along short side 
dimensionless coordinate, y/a 

eigenfunction for case of internal heat generation and qg = q s 

coordinate measured along long side 

dimensionless coordinate, z/a 

dimensionless coordinate, z/b = z’/a 

wall heat flux ratio, q s /qp 

wall heat flux ratio, q s /qg = l/a 

> .> 

eigenvalues for case of internal heat generation and Ip = Is = 
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^tn 

5 n 

0(x',y') 


K 



P 


a 

'Hx 1 , z f ) 


eigenvalues for internal heat generation and q^ = q g = 0 

eigenvalues for vail heat transfer and Q = 0 

eigenvalues for vail heat transfer and Q = 0 

function defined hy equation (32) 

thermal conductivity of fluid 

eigenfunction, 

fluid viscosity 

fluid density 

aspect ratio, h/a 

function defined hy equation (33) 


B 

h 

d 

e 

Q 

q 

s 

v 


Subscripts 

refers to broad vail 
bulk mean value 
fully developed 
entrance value 

insulated vail, internal heat generation 
vail heat flux, no internal generation 
refers to side vail 
value at vail 


* 


entrance region 


Superscript 
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Figure L - Coordinate system for rectangular channel 
with different uniform heating on each pair of op- 
posite walls. 



(a) Uniform heating of broad walls and uniform 
cooling of side walls (a - -1). 


Figure 2. - Wall temperature development In 
thermal entrance region of rectangular ducts. 
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(b) Uniform heat transfer from broad walls only 
(a » 0). Duct aspect ratio, 1 <a<°°. 


Figure 2, - Continued. 



(c) Uniform heating of four walls (a * 1). 
Figure 2. - Continued. 
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Dimensionless longitudinal distance, <x/a)/RePr 


id) Uniform heating of short walls {0 - 0). 
Figure 2. - Concluded. 



(a) Uniform heating of broad walls and uniform 
cooling of sidewalls (a - -1). 

Figure 3. - Longitudinal variation of bulk mean 
temperature for various values of heat flux 
ratio R. 





(b) Uniform heating of broad walls (a - 0). Duct 
aspect ratio, l<o<°°. 


Figure 3. - Continued. 



Dimensionless longitudinal distance, (x/a)/RePr 
(c) Uniform heating of four walls <a - 1). 
Figure 3. - Continued. 
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(d) Uniform heating of sidewalls (n - 0). 
Figure 3. - Concluded 



/ * y/a z" - z/b 

(a) Uniform heating of broad walls and uniform 
cooling of side walls (a = -1). 

Figure 4. - Fully developed wall temperatures of 
rectangular channels. 
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(d) Uniform heating of side walls only (B - 0). 
Figure 4. - Concluded. 



(a) Along dud corner in thermal entrance re- 
gion with uniform heating of broad walls and 
uniform cooling of side walls (a * -1). 


Figure 5. - Wall temperature ratios. 
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Dimensionless longitudinal distance, (x/a)/ RePr 

(b) Along wall centerlines In thermal entrance 
region with uniform heating of broad walls 
and uniform cooling of sidewalls (a • -1). 

Figures. - Continued. 



0 .1 . 2 .3 .4 
Dimensionless longitudinal distance, (x/a)/RePr 


(c) Along duct corner In thermal entrance re- 
gion with heat transferred from short walls 

ra-o). 


Figure 5. - Continued. 
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(d) Along wall centerlines in thermal entrance 
region with heat transferred from short walls 
(0 - 0 ). 


Figure 5, - Concluded. 


